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Background 

Collaborative  effort  \A/ithin  AFRL 

-  Computational  Mechanics  Branch  (RWWC) 

-  High  Explosives  Research  and  Development  Branch  (RWME) 

-  Damage  Mechanisms  Branch  (RWMW) 

Goal:  Predict  survivability  of  energetic  payload  of  high  speed 
penetrating  weapons 

Approach: 

-  UseXCMT 

•  Investigate  internal  structure  of  pristine  and  damaged  materials 

•  Generate  representative  computational  meshes 

-  Perform  meso  scale  simulations  of  composite  energetic  materials  used 
in  AF  penetrating  weapons 

-  Incorporate  deformation  damage  into  material  models  for  energetic 
materials 

-  Extract  data  from  FEM  simulations  of  RVE's  to  feed  thermal-reactive 
prediction  code 
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Status 


•  Established  capability  to  generate  computational 
mechanical  RVE  from  XCMT  data 

•  Generated  experimental  data  for  the  mechanical 
behavior  of  simple  explosive 

—  Pristine 

-  Damaged 

•  Performed  simulations  on  mechanical  RVE's 

-  From  XCMT 

-  Idealized 

•  Developed  and  performed  initial  thermal-reactive 
predictions  (currently  under  development) 
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RVE  Study 
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Pressure  Loads  from 
Systems  Level  Simulation 
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RVE  from  XCMT 
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True  stress  (MPa) 


Characterization  Data  of  Composite 
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True  stress  (MPa) 


Characterization  Data  of  Composite 


Specimen  1  to  5 


-0.1  0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7 

True  strain  (mm/mm) 


Compression  Tests  at  24  mm/min 


Cleared  for  Public  Release  -  96ABW-2012-0263 


Characterization  Data  of  Composite 
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Thermal-reactive  Model 


•  We  can  use  nearest-neighbor  distributions  to 
help  quantify  explosive  morphology. 

•  We  may  estimate  statistical  distributions  of  the 
heterogeneous  temperature  rise  that  results 
from  mechanical  loading. 

•  These  distributions  can  be  used  to  initialize  a 
thermal-reactive  representative  volume  element 
(RVE),  which  is  called  a  "realization". 

•  Many  realizations  can  be  initiated  using  the 
same  statistics. 

•  Not  all  realizations  will  undergo  complete 
chemical  reaction,  even  with  adiabatic  RVE 
boundary  conditions. 

•  We  are  interested  in  this  uncertain  "region",  or 
ignition  threshold.  It  is  our  goal  to  understand 
and  quantify  this  uncertainty. 
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Main  Steps:  Processing  and 
Post-processing 

Step  #1:  Run  and  post-process  the  finite  element  (FE)  RVE 

a.  Specify  a  loading  condition  and  run  the  FE  simulation. 

b.  Post  process  the  following  for  each  crystal: 

i.  centroids,  X, 

ii.  size  or  volume,  V| 

iii.  maximum  temperature, 

iv.  average  temperature,  Tg^g  | 

Step  #2:  Obtain  statistics  (Gaussian  mean  iJ.  and  std.  deviation  a) 

a.  Use  the  cutoff  temperature,  to  classify  each  j  as  either  a  H.S.  (hot  spot)  or  a  C.S. 
(cold  spot). 

b.  For  H.S.  -  T^axj  /  Lvgj/  ^hsj/  (nearest  neighbor) ,  find  ^  and  a. 

c.  For  C.S.  -  Tgyg  j,  j,  and  ,  find  and  a. 

d.  For  H.S.  -  ,  find  find  ^  and  a  for  3  bins  using  all  H.S.  T^gx  values. 
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Main  Steps:  Processing  and 
Post-processing 

Step  #3:  Initialize  and  run  a  realization  of  a  thermal-reactive  RVE 

a.  For  H.S.,  randomly  initialize  centroid  values,  with  VpVhsijU  ,a). 

b.  For  H.S.,  synthetically  anneal  j  to  fit  desired  distribution. 

c.  For  H.S.,  randomly  assign  hot  spot  T^axj- 

d.  For  H.S.,  anneal  the  3  temperature  nearest  neighbor  distributions  using  values. 

e.  For  C.S.,  randomly  initialize  centroid  values,  X^.^  j  with  ,a). 

f.  For  C.S.,  synthetically  anneal  X^.^  j  to  fit  desired  NN^.^  distribution. 

g.  Run  the  thermal-reactive  RVE  simulation. 

Step  #4:  Obtain  ignition  distribution  for  a  given  load,  L 

a.  Repeat  Step  #3  many  times  and  record  ignition  results  (either  yes  or  no  ignition). 

b.  Calculate  the  ignition  probability  based  on  multiple  realizations  for  a  single  load,  L. 

Step  #5:  Obtain  ignition  threshold  as  a  function  of  load,  L 

a.  Repeat  Steps  #1-4. 

b.  Calculate  the  ignition  threshold  as  P|(L)  as  a  function  of  the  load,  L,  where  P,  is  the 
probability  of  ignition  that  is  obtained  from  the  ignition  probability  in  Step  #4. 
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Synthetic  Annealing 


-  Assume  that  the  nearest  neighbor,  NNl,  distribution  is  Gaussian 

P[4||(di  -  5d/2)  <  4  <  (di  +  Sd/2)]  =  ajvjvexp[-i(^)'} 

-  for  i  =  1, ....  Nd,  where  =  number  of  bins  to  discretize  NNl 

-  5d  =  width  of  the  bins 

-  iJ.  =  mean  distance 

-  a  =  standard  deviation  of  the  distance 

-  =  scaling  parameter  so  that  sum  of  di  probabilities  sum  to  unity 


CLnN 
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Synthetic  Annealing 


We  let  P(di)syn  denote  the  synthetic  probability  obtained  from  the 
annealed  configuration  and  P(di)act  denote  the  actual  or  experimental 
distribution. 

The  objective  function  to  be  minimized  is 


^NN  — 


Nd 


2 


Procedure: 

a)  Randomly  assign  crystal  centroids  to  locations 

b)  A  subset  of  the  centroids  are  randomly  perturbed  in  the  x,  y,  and  z 
directions. 
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Synthetic  Annealing 


-  Procedure  (cont.): 

c)  The  objective  function  is  evaluated  again  to  give  a  new 

d)  If  the  new  value  of  's  smaller  than  the  previous  value,  then  the 
new  configuration  is  kept.  If  not,  then  the  probability  that  the 
configuration  will  be  kept  is  given  by, 

P  =  exp{-c1e^/T^.} 

Ty  =  virtual  temperature  (can  be  a  function  of  number  of  attempts) 

e)  Steps  b-d  are  repeated  until  either  i)  the  maximum  number  of 
attempts  has  been  reached,  or  ii)  the  number  of  temperature  (TJ 
decrements  has  been  reached. 
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Example:  Step  #3 


The  random  centroid  distribution  is  shown  below. 
Hot  spots  (hs,  red)  and  cold  spots  (cs,  blue). 

=  4  to  12  (varies  linearly  with  load) 

Nt„tai  =  30  crystals 

^cs  “  ^total  “  ^hs 

Lx  =  Ly  =  Lz  =  15  microns 
/J(D,hs)  =  M(D,cs)=  2  microns 
'^(D,hs)='r,D,cs)=l  microns 
Min  diameter  =  0.5  micron  (hs  and  cs) 

M(NNi,hs)=  (4-5  -  Nh3/20)*le-6  microns 
i^(NNi,cs)=  (3-0-  N,3/20)*le-6  microns 
^(NNi,hs)“  ^(NNi,cs)“  ^  microns 
nx  =  ny  =  nz  =  31 


15  0 


y, 


X,  fim 
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robability 


Example:  Step  #3 

Hot  spot  centroid  annealing  -  nearest  neighbor  distribution(NNl) 

^attempt  “  300,  -  2*dx,  Tyj^ual  “  ^move  ~  0-20 


Hot  3pot3:  initial  NN1  Hot  3pot3:  anneaied  NNl 


Distance  between  centroids,  /um  Distance  between  centroids,  /um 

ENNl  (random)  =  0.4151  ENNl  (annealed)  =  0.1679 
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Example:  Step  #3 

Non-hot  spot  centroid  annealing  -  nearest  neighbor  distribution(NNl) 


Cold  3pot$:  initial  NN1 


Distance  between  centroids,  (im 


012345678 
Distance  between  centroids,  /tim 


ENNl  (random)  =  0.3837 


ENNl  (annealed)  =  0.0939 
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Example:  Step  #3 


Demonstration  assumptions  for  this  example: 

-  All  grains  have  the  same  volume. 

-  Initial  temperature,  within  each  grain,  is  equal  to  the  bulk  temperature  of 
the  RVE. 

-  The  material  properties  Cp,  rho,  kc,  and  dHrxn  are  all  constant. 

-  The  hot  spot  maximum  temperature  distribution  is  constant. 

-  Loading  is  simulated  by  linearly  varying  the  bulk  temperature  and  the 
number  of  hot  spots. 

-  There  is  no  correlations  between  crystal  temperature  and  the  NNl 
distribution. 

-  The  RVE  is  adiabatic,  and  it  has  symmetric  boundary  conditions. 
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Example:  Step  #3 


(Right)  Temperature  contour  on  one  face  of  the  RVE  cube  at  300  microseconds. 
Periodic  boundary  conditions. 

Run  time  =  300  microseconds 


M(xhs)(avg)  =  800  K 
M(T,buik)(avg)  =  522  to  625  K 
^5'(T,hs)(avg)  =  100  K 
^^(xcs)(avg)  =  100  K 
CFL  =  0.6 

kc  (conductivity)  =  0.2  W/m*K 

Cp  =  1466  J/kg*K 

Rho  =  1180  kg/m^3 

dHrxn  (heat  of  reaction)  =  2.0e9  J/m'^3 

Tcriticai  =  ^  (crystal  ignition  criterion 


X,  microns 
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Example:  Steps  #4-5 

A  single  loading  condition  shows  the  average  crystal  temperature,  not  the  bulk  RVE 
temperature. 

20  realizations  for  each  loading  condition. 

Threshold  is  approximately  530  to  580  K  (50  degree  range) 


Average  crystal  temperature  Ignition  threshold 


Cleared  for  Public  Release  -  96ABW-2012-0263 


Summary 


Can  generate  meshes  from  XCMT 

Have  experimental  data  of  mock  and  explosive 
composite 

Have  carried  out  simulations  on  mechanical 
RVE's  from  XCMT  and  idealized 

Developed  and  demonstrated  thermal  ignition 
predictive  approach 
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Future  Work 


•  Develop  approaches  for  extracting  data  from 
FEM  to  feed  thermal  ignition  code 

•  Improve  temperature  computation  approach 
in  FEM 

•  Implement  better  thermal-reactive  modeling 

•  Quantify  the  range  of  validity  using 
experiments 
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Related  Efforts 


•  SBIR's 

-  Phase  2  PENETRATION  SURVIVABLE  ADVANCED  ENERGETIC 
—  Symmetrix 
—  Streamline 
—  Illinois  RocStar 

•  STTR's 

—  Phase  I  with  Multiscale  Design 
—  Topic  on  NDE  for  low  density  gradient  materials 

•  Los  Alamos  NL 

•  Should  soon  have  a  signed  Project  Agreement  with 
Germany 

•  Udaykumar  of  University  of  Iowa 
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